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Abstract 

Recent results from the PAMELA experiment indicate an excess in the positron spectrum above 
10 GeV, but anti-proton data are consistent with the expected astrophysical backgrounds. We 
propose a scenario that reproduces these features. Dark matter annihilates through channels 
involving a new heavy vectorlike lepton which then decays by mixing with Standard Model leptons. 
If charged, this heavy lepton might be produced at the LHC, and could lead to multi-lepton final 
states or to long-lived charged tracks. Large neutrino detectors such as ANTARES or IceCube 
might be sensitive to a monochromatic neutrino line. This scenario may be simply embedded in 
various models, including an extension to the NMSSM. 
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I. INTRODUCTION 



There is great deal of excitement surrounding recent measurements of the positron flux 



by the PAMELA experiment 



2] , which show an increase in the positron fraction above 



10 GeV. Interactions of cosmic ray protons with the interstellar medium would be unable to 
produce a significant rise, prompt 
well as a reexamination of pulsars 
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5[, in light of earlier predictions f|. 
Interpreted as a dark matter signal, the PAMELA result places strong constraints on 
annihilation channels. The sharp rise indicates a hard lepton spectrum [?], Q , while the lack 
of an accompanying rise in the p/p ratio strongly constrains hadronic annihilation modes 
0. Consequently, if PAMELA results are taken at face value, the dark matter particle 
annihilates primarily to leptons — hadronic channels are suppressed. 

Construction of a dark matter candidate that primarily annihilates directly to leptons is 
non-trivial. Majorana fermions have helicity suppressed annihilations to light fermions, and 
Dirac fermions face significant model building challenges 10[] . However, one can consider 
new annihilation modes into non-Standard Model states, which then subsequently decay 
to leptons. Such models can be divided into two categories. In the first approach, dark 
matter annihilates to light intermediate particles. If sufficiently light J~ GeV) , the absence 
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of anti-protons can be explained by simple kinematics 11[ (see also 12 
second approach is for the dark matter to preferentially couple to leptons via dynamics jl7l |. 
Alternative models have considered the pos sibility where dark matter decays preferentially 



to leptons 18, 
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Here we present a hybrid approach. Dark matter dominantly annihilates to final states 
involving new, heavy particles. These particles mix with the Standard Model leptons, and so 
carry lepton number. Following the decay of these new particles, dark matter annihilations 
produce hard leptons in association with gauge bosons. The annihilation spectra from this 
model are distinct from those previously included in the literature thus far. 

In the next section we discuss the general framework for dark matter annihilations and 
heavy lepton decays. The scenario is a very simple extension of the Standard Model. In 
section III, we calculate the cosmic ray fluxes from dark matter annihilations and the boost 
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factor relative to the thermal cross section. We find this scenario can explain the PAMELA 
excess with a small boost factor. In section IV, we discuss the possible observation of 
high energy neutrinos from dark matter annihilations. In section V, we discuss the possible 
production of new vector leptons at the LHC, and their signatures. In section VI we construct 
extensions of the Standard Model and Supersymmetric Standard Model which realize this 
scenario. Finally, we conclude. 



II. ANNIHILATIONS OF DARK MATTER INTO A NEW, HEAVY LEPTON 

The PAMELA data do more than confirm the excess seen by HEAT [2^] and AMS-01 



281 ] . The sharp upturn in the positron fraction seems to disfavor annihilations with hadronic 
modes, because the positrons from a hadronic shower are generally too soft to produce such 
a spectrum. Furthermore, such showers contain anti-protons, and no excess is observed 
in this channel. Thus, to reproduce the positron/electron ratio and the antiproton/proton 
ratio observed by PAMELA, dark matter annihilations should produce a larger ratio of high 
energy positrons to antiprotons than annihilations to vector bosons [63]. 

We propose a new mechanism to achieve this. We assume dark matter annihilation 
proceeds via the processes in Figure [TJ The dark matter particle, Xi annihilates through new 
channels which involve heavy vectorlike states, H/5, which carry lepton number and thus 
subsequently decay into light leptons and W or Z bosons. 

This theory can explain the PAMELA signal. The introduction of new, heavy leptons 
allows direct annihilation into leptonic states which are not suppressed by a small fermion 
mass. This can occur either from direct annihilations to heavy leptons alone (i.e. through a 
XxHS-type operator) or from a combination of new heavy leptons and light leptons (through 
a xx^rtype operator). Because we do not rely upon a helicity flip on the external line, 
such processes can be large and dominate the annihilation rate. We will not specify the 
mechanism by which the annihilation proceeds for now, but will revisit this question. (A 
similar annihilation to a heavy/light pair can occur in RS theories, although with different 
masses and spectra than we consider here |29|.) 

Under SU{2) x U(l) the S can have quantum numbers (2,1/2) (like a Standard Model 
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FIG. 1: Left: The annihilation process of the dark matter x to a heavy S lepton and either S or a 
Standard Model lepton I. The shaded circle represents a model-dependent annihilation mechanism. 
We discuss possibilities for this in section IVI1 Right: Decay of the H into to a lepton and a W 01 
Z boson. 

£), (1 , 1) (like a Standard Model e c ), or (1 , 0) (a complete Standard Model singlet, n). Any 
of these is viable, but they have distinct phenomenology. 

In the case where we add vectorlike particles with quantum numbers of e c or n, decays 
proceed with comparable rates both through W and Z bosons [64]] . Thus, for XX — — ■ each 
annihilation (on average) roughly yields a single charged lepton and two gauge bosons. In 
contrast, as we shall see, annihilation to a pair of vectorlike heavy doublets ultimately yields 
final states with charged leptons and gauge bosons in a one to one ratio. Since gauge bosons 
are a source of anti-protons, taking H as a vectorlike SU(2)l doublet produces the highest 
e + jp signal, as is desirable to fit the PAMELA data. 

For the study of cosmic ray signals, then, we will focus on the £-like case. For the purpose 
of the PAMELA data, the consequence of considering the n- or e c -like models instead is to 
multiply the p spectrum of the case at hand by a factor of ~ 2, as an increased boost factor 
would be needed to fit the positron data (see section UTT|) . 
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For concreteness, we begin with a Lagrangian|65] 



(1) 



Prior to electroweak symmetry breaking, we can do a rotation in S) space to remove the 
€i terms at the cost of modifying the & terms. Subsequent rotations on the £j and el allow 
diagonalization of the 3x3 submatrix. After electroweak symmetry breaking and performing 
the above rotations, the charged and neutral fermion mass matrices are 

/ m fi \ / \ 





M 



charged 



m M 
m T 



M 



neutral 





W 



(2) 



Physical neutrino masses are sufficiently small that we can neglect them. There is no a priori 
reason why m e < m^^rriT. In anarchical modelsfs^], they may all be the same, or should H 
carry a lepton flavor charge, it could even be that m e ^> m M , fh T . For future discussion, it is 
useful to introduce the small parameter Si = fhij ' \x. 



A. Decays of the heavy leptons 

Let us begin by focusing on the phenomenology of the neutral sector. All neutral states 
(both Standard Model and the new heavy state) have identical couplings to the gauge bosons. 
So, when the mass matrix is diagonalized, no off-diagonal couplings to the Z boson are 
created — the mathematics is identical to that which ensures the absence of FCNCs in the 
Standard Model. However, the differences in diagonalizing the charged lepton mass matrices 
will lead to off-diagonal charged current decays. Thus, neutral E states decay to a W ± t^ 
final state. 

The decay modes of S 1 * 1 are more complicated. We assume that the S 1 * 1 state is heavier 
than the S° state. The charged state has non-vanishing couplings to three final states: 
S°iy ± * , Z^, and W ± v. Typically, the most suppressed of these decays is the decay to 
W ± v. This heavy-light charged current coupling is down by a factor (m^//i)5j, where is 
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the relevant charged lepton mass. Were to vanish, then the rotation on the left-handed 
doublet has no effect (just as in the neutrino case above). This is the origin of the suppression. 
Thus, the two decay modes of potential phenomenological interest are E W* ± and Z£ . The 
first of these receives no suppression from small couplings, but receives a large phase space 
suppression. If no additional couplings are added to the theory, the dominant contribution 
to the Am= = m s ± — m^o mass splitting is from the loop effect induced by the Coulomb 
interaction, analogous to — 7r° mass splitting. Because this effect is cut-off by the scale 



of electroweak symmetry breaking, it is finite and is approximately given by Am= 



aM z 
2 



350 MeV. There is an additional (subdominant) contribution to the splitting between the 
charged and neutral states of size 5 2 fi due to the mixing with the Standard Model states. 
We will see later that S is constrained to be < 10~ 2 . For the tiny splittings induced by the 
Coulomb interactions, decays will proceed to 2° and a soft charged pion with a rate given 



by 



3l| 



I\ 5 o = G|cos 2 ^^ 2 (Am s ) 3 t /l - f^-j . (3) 

The remaining decay, E* 1 — > Z£^, is unsuppressed kinematically, but the relevant coupling 
is proportional to 5. It is given by 

m~<5 2 g 2 



647r cos 2 9 



w 




(4) 



For 5 > 1CT 7 , the kinematic suppression of the charged current decay is more severe, and 
the decay via H 1 * 1 — > Zi^ dominates. In the remainder of the paper, we will assume that 5 
is sufficiently large so that S 1 * 1 decays in this way. 



B. Spectrum of annihilation products 

We first consider the case when the dark matter dominantly annihilates into two heavy 
states. It is natural to assume the dark matter annihilates equally into and 
Based on the above discussion, these processes yield £ + £~ZZ and £ + £~WW final states. In 
principle, the mixing of the new vectorlike state can be with any Standard Model lepton 
family. Not surprisingly, the case where mixing is dominantly into electrons provides the 
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most favorable to fit the PAMELA positron spectrum: it produces the hardest positron 
spectrum possible for a certain amount of antiprotons. Mixing with \i and r leptons will 
require larger boost factors to produce the positron signal and will result in some tension 
with the anti-proton data. 

We now address the spectrum that results from dark matter annihilation. It depends 
on the masses m x and mn. In the limit m= — > m x , the spectrum decomposes into two 
components: gauge bosons of energy ~ m x /2 (assuming m= ~ m x 3> my) and electrons 
and positrons of energy m x /2. Here my is the mass of the relevant gauge boson. In the 
opposite limit, m x —>■ my, the energy of the positron from the 5 decay goes to zero. The 
most relevant positrons come directly from gauge boson decay will have maximum energy 
m x . Finally, one finds the hardest electron/positron spectrum for the lepton from S decay 
comes for 



The spectra for these three limits are shown in Fig. [2j 

The addition of two hard leptons to every annihilation allows a fit to the observed positron 
fraction with a smaller boost factor than would be required if annihilation proceeded directly 
to gauge bosons. This provides a relative suppression of the antiproton contribution. This is 
roughly an order of magnitude change with respect to the case of pure annihilation to gauge 
bosons. Residual p's will still provide a constraint. We will explore this in section HTT1 

Next, we consider the case where dark matter annihilates into one heav y a nd one light 
lepton. Similar final states occur in Randall-Sundrum theories of dark matter [29] . Again, the 
natural assumption is that dark matter annihilation populates the charged and neutral states 
equally. As long as m= < 2m x , the H final state is available. In this limit (mn — > 2m x ), no 
kinetic energy will be available for the S£ final state. Thus, the monochromatic light lepton 
is soft. Effectively, the final state is a S field at rest, which then decays as described in the 
previous section. In the opposite limit of a very light S, there is a combination of a hard 
monochromatic lepton plus a H decay spectrum. We plot the positron injection spectrum in 
three cases, win = 95 GeV, ~ m x , and ~ 2m x in Fig. [3] for the case where the final state 
lepton is an electron. In this heavy-light case, a monochromatic neutrino is also produced. 
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FIG. 2: Positron spectra from the process XX ~^ for three different E masses for a dark matter 
mass m x = 800 GeV. These injection spectra are shown prior to propagation and correspond to a 
single annihilation. It is assumed that the E decay via a mixing with electrons. 

We will discuss its consequences in Section IV. 



C. Precision Constraints 



Although we introduce new states that mix with Standard Model leptons, it is obvious 
from the outset that there need be no conflict with precision observables. If the mixing 
with the Standard Model (as controlled by 5 = rh/fi) is sufficiently small, there will be 
no measurable effect on, e.g., g — 2, although S will still decay on (cosmo logically) short 
timescales. However, to understand the phenomenology of the heavy states, it is important 
to consider how large 5 can be. 

In the case of a vectorlike doublet, there are two relevant classes of constraints. First, if 
both 5^ and 5 e are non-zero, a bound arises from the process fi — > 3e. Consistency with the 
non-observation of this process BR(n — > 3e) < 10~ 12 [32j forces ^S e 5^ < 2 x 1CT 3 . A weaker 
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FIG. 3: Positron spectra from the processes XX ~* ^ e f° r a dark matter mass m x = 800 GeV and 
three different H masses. These injection spectra are shown prior to propagation, and correspond 
to a single annhilation. 

constraint arises on this combination from fi — > ej. Second, there are constraints present 
even if only one S is non-zero. In particular, there are precision electroweak constraints from 
considering, e.g., the universality of T(Z —>■ It). This forces 8g < 2 x 10~ 2 

Should we choose the vectorlike state to instead have the quantum numbers of a right- 
handed electron, there is also an important bound that comes from the universality of Gf- In 
that model, the SU(2) ^-singlet exotic is mixed with a standard model lepton and effectively 
reduces the coupling of this light state to the W. Thus, for the singlet model there is an 
additional constraint, 5e < 4 x 10~ 3 . 

So, there is a large window (roughly 10 -7 < 5 < 10 -2 ), where the decay phenomenology 
described in the previous section applies and is consistent with all precision constraints. We 
assume that we are in this window and now move to a detailed discussion of the astrophysical 
signals of this model. 
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III. EXPLAINING THE PAMELA RESULTS 



While we have presented a simple scheme for generating hard positrons with few 
antiprotons, it remains to be seen whether we can fit all of the existing data, and what 
the minimal boost factor above the thermal cross section is required to reproduce a good fit. 

To investigate the annihilation spectra, we implemented the processes described in section 
HT1 into MadGraph [33( . This allowed the simulation of the annihilation and decays shown 
in Fig. [TJ Subsequent decays of the W and Z bosons and hadronization where done with 



PYTHIA 



34|. This produced an injection spectrum of positrons, electrons, antiprotons, and 



gamma rays. This spectrum was input to GALPROP 



351 ] to pr opag ate the decay products 



36) 



through the galaxy. For concreteness, an NFW profile was used 

with r s = 20 kpc, where r Q = 8.5 kpc is the galactocentric distance of the sun and p Q = 0.3 
GeV/cm 3 is the local dark matter density, although the shape of the profile will have little 
effect on the observed positron fraction. For positron, electron, proton, and antiproton flux 



backgrounds, we use the model described in 11] with the Alfven velocity va = 20 km/s. The 



diffusion coefficient is taken to be D = (3(5.88 x 10 28 cm 2 /s) (R/4 GV) 1 / 3 , where f3 = v/c 
and R is the rigidity. The height of halo region is set to L = 4 kpc. 

There are a number uncertainties that enter into the calculation of the rates. In addition 
to the cross section itself, indirect detection signals are proportional to the number density 
squared of dark matter particles. The uncertainties associated with this are usually encoded 
into a "boost factor" . We define the boost factor as 

BF=J-([ A '^M < OT > ) . (7) 

V C R \J n 2 (r) {av} therma i J 

Here n(r) = p(r)/m x is the previously mentioned NFW profile with local density po — 
0.3 GeV cm -3 . n true {r) = ptrueij)jm x is the actual (possibly clumpy) number density of the 
dark matter particles in the halo, and the integration is over a region Vcr defined by cosmic 
ray propagation. 

The integral over the number density contains all of the information on astrophysics 
and is subject to both uncertainty in the halo profile and uncertainties on propagation. 
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dumpiness gives an increase in n 2 when compared to a smooth profile. Should the boost 
factor arise from astrophysics, it is expected that the boost factors could be different for 
positrons and antiprotons. Depending on the positions of subhalos with respect to the solar 
position, it could lead BF e +/BF p ~ 3 and energy dependent boost factors. We will 
ignore the possibility of energy dependent boost factors in this work, but will bear in mind 
the possibility that astrophysical uncertainties might somewhat ameliorate tensions between 
positrons and anti-protons. Note that uncertainties in the local number density could be a 
factor of ~ 2 (yielding a boost of ~ 4) [38|, |39j. Altogether, boosts of up to ~ 10 arising 
from astrophysics alone are likely reasonable, while much higher boosts (as we shall find are 
necessary for heavier WIMPs) likely would rely on a contribution from a non-trivial dark 
matter cosmology. 

It should also be noted that uncertainties in the propagation of antiprotons could suppress 
the dark matter contribution to the antiproton signal by up to an order of magnitude 9|. 
The effects of these changes on the positron flux is more modest. So, while in this work we 
quote a boost factor for the positron signal and naively apply the same the boost factors to 
the anti-protons, it should be kept in mind that masses that initially appear to be ruled out 
by the observed p/p ratio might be allowed once a full accounting of these uncertainties is 
made. 

We show the results of our analysis in Fig. |H We find that it is possible to fit the 
PAMELA positron data for a variety of masses. For each set of masses, we find the boost 
factor (as defined above) that gives the best fit to the data. To avoid the complications 
of solar modulation {40! . I4I ]. we perform a x 2 fit to only the four highest energy bins of 
the PAMELA data (where solar modulation is known to have little effect). The results 
are the boost factors shown in Fig. HI Since PAMELA data were taken during a negative 
polarity part of the solar cycle, correcting for modulation would reduce positron ratio at 
lower energies, bringing the curves into qualitatively better agreement with the data. 

The smallest dark matter mass found to be consistent with the PAMELA positron data 
is m x = 200 GeV with m= ^> mw- In this case, the boost factor required is only 6, which 
could very plausibly be entirely due to astrophysics. Thus, dark matter signals could result 
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from a thermal annihilation cross section with no need for non-thermal production or late 
time cross section enhancements. As m x increases, the boost factor required is increasingly 
unlikely to be given by astrophysics alone, and an increase in the annihilation cross section 
beyond the thermal one is likely necessary. This could be consistent, for example, with a 
cosmological history where the dark matter is produced non-thermally, or if the late-time 



properties of the dark matter change 



42( | . If H decays to muons or tau leptons, this will 



produce a softer positron spectrum and will require boost factors to fit the positron data. 
These will come into more tension with the antiproton data. For muons, the difference is 
roughly a factor of two from the electron case. 

Dark matter masses m x ~ 500 GeV are in tension with the antiproton data. Keeping 
in mind the earlier caveats mentioned on astrophysical uncertainties (both propagation and 
dumpiness), a prediction of this model is that the dark matter have mass less than 500 GeV. 
A turn-over or plateau in the positron fraction should be imminent. A reasonable fit to both 
the p and positron data is possible in the case where the 5 dominantly mixes to muons for 
m x < 400 GeV. 

When annihilation occurs via XX ~^ the positron spectra have slightly cuspier shapes 
than when annihilation proceeds through XX ~^ EE, see Fig- El The boosts required are 
similar to the EE case, and the fits are comparably good. Annihilations via XX ~^ E/i can 
also give good fits, but require somewhat larger boosts, leading to some tension with anti- 
proton results, even for lower dark matter masses. A dark matter mass of m x = 200 GeV that 
decays in this way is marginally consistent with the positron and antiproton observations 
for a boost factor of 10-20, depending on the S mass. 

Finally, if XX ~~ EE annihilation dominates, a choice of m x = 800 GeV and m= = 370 
GeV allows a qualitative fit to the shape and position of the ATIC data 43J], demonstrated 
in Fig. [6j In this figure, the background is scaled such that it saturates the low energy data. 
While the signal does not reproduce the ATIC data precisely, we find this qualitative fit 
encouraging. If the ATIC anomaly persists, and this model is to explain it, this points to a 
m x = 800 GeV to 1 TeV dark matter particle with a large boost factor, ~70. Non-standard 
choices for the propagation parameters would be necessary in to reconcile the high energy 
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FIG. 4: Positron fraction and antiproton fractions XX ~^ EH. The boost factors are relative to 
3 x 10 26 cm 3 /s, and are found by scaling the positron flux to fit the observed positron fraction. In 
the m x =200, 300 GeV cases, the boost factor required is sufficiently small that the dark matter 
might have a thermal relic abundance. 
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FIG. 5: Positron fraction and antiproton fractions for \X ~* The boost factors are relative to 
3 x 10 26 cm 3 /s, and are found by scaling the positron flux to ft the observed positron fraction. In 
the m x = 200 GeV case, the boost factor required is sufficiently small that the dark matter might 
have a thermal relic abundance. 
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FIG. 6: ATIC signal for m x = 800 GeV and m~ = 371 GeV. The background has been scaled 
to saturate the low energy total flux. The required boost factor is 72, indicating a non-standard 
cosmological history for the dark matter. 

antiprotons production with the PAMELA data (Fig. Hj). 
IV. NEUTRINO SIGNALS 



The detectability of neutrino signals, in particular as related to models to explain 



'AMELA has recently been studied 



44 



45 



461 ] . Upcoming experiments such as IceCube 



^ and ANTARES 48] could potentially show sensitivity. 

If annihilation proceeds as XX ~ * EE, with E an electroweak doublet, there are essentially 
no hard neutrinos produced. However, in the case that annihilations proceed XX ~^ ^2, such 
processes will produce equal numbers of monochromatic charged leptons and neutrinos. This 



monochromatic v presents an exciting experimental 



The background from cosmic ray neutrinos 49] can be well modeled by a power law. 



target. 



Following 44j, we examine the signal at a neutrino telescope by requiring a signal to 



background ratio of 1 in a bin from 10 5 m x to m x when compared with the angular average 
of Vn, + v a . We assume that oscillations yield 1/3 of all neutrinos as muon flavor. The signal 
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is proportional to the line-of-sight integral of the dark matter density averaged over the 
relevant solid angle, Jah- We calculate Jac of 680 for a 2° region around the galactic center 
in an NFW profile. This solid angle corresponds to the angular resolution of ANTARES. 
Taking this gross analysis requiring S:B=1, we find that a detectable signal can be found for 
boosts of roughly 320 66]. 

Such a limit depends only weakly on the dark matter mass (assuming fi <C rn x ), although 
for very large neutrino telescopes where the background statistics are large enough to model 
the background, the appropriate comparison is naturally S/ VB and greater sensitivity would 
be possible. For reference, in ten years of IceCube running we expect 200 background events 
from 500-1000 GeV in a 2° region [46]. Thus, without a careful examination of background 
models, we cannot say how detectable this scenario would be at lower masses. For higher 
masses, detection will be statistics limited, and the S:B criterion is the right one. This case 
seems borderline in its detectability. Ruling out this model with neutrinos seems difficult. 
That said, there is substantial evidence within numerical simulations for many subhalos 



in the Milky Way[50j, |5l|, |52j. Should one of these provide a larger J An, for instance by 



being cuspy or very nearby, a positive detection remains tantalizingly possible. IceCube has 
good angular resolution and could focus on one of these regions in an effort to suppress the 
atmospheric neutrino background. 

In fact, if this substructure is sufficiently generous, then it might be possible to observe 
the z/'s produced in the case where the S is an SU(2)l singlet, but dark matter annihilations 
go to HH. In this case, the neutrinos have a spectrum very similar to that seen for positrons 
in Fig. [2j They have energy roughly a factor of two lower than the case where direct Sz/ 
dominates. Due to the steeply falling power law nature of the background, these neutrinos 
would have to be visible over a background roughly an order of magnitude higher. 

Lastly, we should note that the observability of the galactic center signals rely on cuspy 
profiles into the inner 2°. Should those profiles extend into the inner ~ 0.5J\ strong gamma 



ray limits from HESS may constrain possible neutrino signals 53|, |54j, |55|, l56j ]. 
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V. COLLIDER SIGNALS 



While dark matter signals are often elusive at a collider, we have the exciting possibility of 
directly producing the new H states that dominate dark matter annihilation. The cases where 
the S has quantum numbers of a leptonic doublet or a right-handed electron are interesting. 
These states will be produced via the Drell-Yan process. The collider phenomenology of the 
case where the H is a complete Standard Model singlet is uninteresting, simply because it 
will not be produced at an important level. 

First, we consider the doublet case. Production of and [E^S are possible. 

The production of leads to a ZZt^l" final state. The cross section at the LHC for 

this mode is 27 fb at m-= = 200 GeV. In 300 fb -1 , we expect ~ 50 six lepton events, with 
4 leptons reconstructing to a pair of Z bosons. While very few events are expected, they 
would be spectacular; irreducible physics backgrounds in this channel are vanishingly small. 
In addition, one could expect observation in the [ZZii — > A£ + 2j) mode and perhaps even 
ihe 2£ + 4j mode. We estimate the background in the 4£ + 2j case to be 9 fb from ALPGEN 



571 ] using the process ZZjj — > 4£jj. These backgrounds could be substantially reduced by 



requiring that the jets reconstruct a Z and applying a Z veto on one of the pairs of leptons. 
Lepton flavor tagging could also additionally suppress backgrounds. Before cuts, we expect 
~ 500 4£+2j signal events in 300 fb~ x for m~ = 200 GeV. Assuming these cuts reduce the 
background by a factor of a few without compromising the signal, this will give a S:B~1. 
The signal would fall to ~ 180 events in 300 fb -1 at m= = 300 GeV. At the Tevatron, the 
cross section for production for m s = 200 GeV is 1 fb, too small to be seen. 

If the mixing between the S ± and the Standard Model is very weak, then the H could 
be long-lived on detector time scales. However, in the doublet case, the existence of the 
S ± — > E°ii ± decay mode unsuppressed by 5 (Eqn. [3]), means that the tracks do not extend 



much more than a centimeter [3JJ. This makes observation of these events unlikely. This 
changes in the case where the E is a SU(2) singlet. Then for sufficiently small mixings, a 
striking long-lived charged track is possible, as all 5 decays are suppressed by 5. Depending 



on the size of S, there is also a possibility 58|, |59j that H could solve the primordial lithium 



problem. For a choice of S that allows S ± to decay around BBN, S ± would clearly be long- 
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lived on detector time scales. The current bound on the production cross section of such a 
long -lived particle from the DO experiment is < 10 fb[60]. This translates to a lower bound 
on the mass of a S'?7(2)-singlet S 1 * 1 of 120 GeV. 

Returning to the doublet case, we turn to Drell-Yan production of This produces a 

W + W~£ + £~ final state. The cross section at the LHC for Drell-Yan production of the 2° is 93 
fb at m a = 200 GeV. In 300 fb -1 , we expect ~ 1300 A£+ # T events. While Standard Model 
diboson production presents an important background, lepton flavor tagging could provide 
a helpful discriminant: the Standard Model rate for 3££' is very much suppressed. It might 
also be possible to observe this final state in the 3£ + 2j+ $t channel. The Standard Model 
background (crx BR) for this state from WZjj is roughly 200 fb. So, it will be important 
to impose a Z veto. Fakes from top production could also be important, and more study is 
warranted. 

Finally, E^H can be produced via an intermediate W boson. This state leads to a 
ZW ± £ + £~ final state, which can ultimately yield striking 5£ events. For m~ = 200 GeV, the 
production cross section is 756 fb. For m= = 300 GeV, the production cross section is 176 
fb, and for m= = 600 GeV, the production cross section is 12 fb. At the Tevatron the cross 
section for [E^S production for m= = 175 GeV is 76 fb, and m= = 200 GeV is 4 fb. For the 
lowest masses, a few multilepton events could be observed. 

VI. MODELS 

One of the appealing aspects of our proposal is that it is simply incorporated into the 
Standard Model and almost any model of new physics. In this section we will outline more 
concrete models that might realize this scenario, first in extensions of the Standard Model 
then in extensions to the minimal supersymmetric Standard Model (MSSM). In each case 
the basic building blocks are the same: namely, the new leptons states themselves, as well 
as an annihilation mechanism, which is accomplished with the introduction of a new scalar. 
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A. Non-supersymmetric implementation 



We begin by considering a non-supersymmetric realization. The vectorlike heavy lepton 
mixes with the Standard Model leptons as in Eq. (JTJ), and decays via a light lepton and Z 
boson. The Lagrangian is 

£ = VxX 2 S + ys^S + ye R ^h + H.c. - V(S), (8) 

where we have rotated away a potential mass term £EH. We also introduce a complex scalar 
S = s + ia with potential 

V(S) = m 2 s S^S + [3S^S 2 + p*SSV + Xs{S^S) 2 , (9) 

which has a global U(l) symmetry explicitly broken by (3. If we allow a vacuum expectation 
value for S to give mass to the dark matter and H, then a is a pseudo-Goldstone boson 
with its mass controlled by the symmetry breaking parameter (3. We have neglected a 
possible S^SH^H mixing term. The presence of this term could affect the pseudoscalar 
phenomenology. While this term will be generated at the loop level, the radiatively generated 
size is suppressed by two powers of rh/fi and a loop-factor - small enough that it has no 
significant effects. 

We consider the case where the dark matter annihilates through the pseudoscalar a . This 
avoids the suppression by the velocity of the dark matter v 2 that occurs in the case of scalar 
mediation. The annihilation cross section is 



2 2/ 2 

Vxm 1 1 m i 



( av )^ = -°x»B-Jl a _ (io) 



In the limit that rn^/m 2 , << 1 and 2m= ~ m x , 

. , ( 94 cm 3 \ , , /300GeV\ 2 

(av) 3E ~ (^5.56 x 10- 24 ^J y\yl J . (11) 



Annihilation into pseudoscalars in this model could also be important, depending on 
the structure of the model. If this occurs, the annihilation into a with subsequent decays 



a — > e + e could provide a good fit to the data [Hi, [12J. However, in the present realization, 
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FIG. 7: Decay of the S into to a lepton and a pseudoscalar a, which subsequently decays to leptons. 

annihilation into light pseudoscalars is suppressed by v 2 since Majorana fermions can only 
annihilate into CP-odd states. 

A simple modification of the above Lagrangian would realize the heavy/light annihilation 
scenario of the previous section. For example, suppose S couples as SG, and SS gets a mass 
from some other source. S no longer gets a vacuum expectaion value, but can still mediate 
the Dark Matter annihilation. 

For light pseudoscalars, a new decay channel S 1 * 1 — > ae^ opens, shown in Fig. [7J The 
pseudoscalar subsequently decays to electrons. The decays of H° to pseudoscalars vanish in 
the limit of vanishing neutrino masses. The branching ratio for charged H decays into light 
pseudoscalars is given by: 



(1-4)2 

v mi. ' 




This branching fraction is small if m a — > m= or if mn >> mz- So for a low mass H, the 
branching fraction to the pseudoscalar will dominate over that to Z bosons until m a ~ mz- 
If this occurs, this will induce a slight modification of the results of the previous section. 
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The net result of S 1 * 1 — > ae dominating is a reduction in both the necessary boost factor and 
the antiproton contribution. Only the S° — > We decays contribute to the antiprotons in this 
case. This light pseudoscalar will have a very small coupling to electrons proportional to b\ 
and thus can easily evade detection at LEP. 

B. A Supersymmetric Case 

The supersymmetric generalization is straightforward. We will assume there is a dark 
matter sector containing the dark matter field X, where the fermionic component is the 
dark matter, that is coupled to the MSSM by a singlet S. This singlet then couples to a 
vectorlike heavy lepton. The superpotential is given by 

W = W M ssm + W DM + W s , (13) 

where the dark matter sector superpotential is 

W DM = VxXXS + y s SEE, . (14) 

In this case, the dark matter annihilates via s-channel exchange of the the singlet field, S to 
a vectorlike lepton. Annihilation via t- and u-channel processes to the singlet S can also be 
important depending on the relatvie sizes of y~ and yx- 

Should S have the quantum numbers of a right-handed neutrinos, we have the mixing 
superpotential terms 

W E = y£EH u + fiEE. (15) 

Then the E marries mostly S with a small amount of of the neutral component of the doublet 
£. The S will decay via mixing with the neutrinos, so possible decays are to the Z boson and 
a light neutrino or to a light lepton and a W boson. Decays of the S into a light neutrino and 
a Z boson will produce the cosmic ray spectra of decays of Z bosons. Decays to leptons will 
produce a hard lepton and the cosmic ray spectra from W boson decays. Essentially identical 
remarks apply to the analogous superpotential where we give S the quantum numbers of the 
right-handed electron. 
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We now focus on the case where the E has the same quantum numbers as the left-handed 
lepton doublet, most similar to the case discussed in the previous section. After appropriate 
rotations, the leptonic portion of the superpotential can be brought to the form. 

W s = y£e R H d + yeEH d + fxEE. (16) 

This superpotential recovers the scenario of the previous section, assuming that the X 
particles are the dominant component of the dark matter, and possess a I2 symmetry in 
addition to the usual R-parity of supersymmetric models. For the X fields to dominate the 
dark matter density, we will need to assume that the neutralino of the MSSM sector is a 
subdominant component. For example, this could be a light thermally produced wino. 

Embedding in a supersymmetric model may lead to possible new signals at the LHC The 
LHC will copiously produce particles with SU(3)c quantum numbers, which will cascade 
down to the MSSM lightest supersymmetric particle. Should cascades contain a neutralino 
with a large bino component, this neutralino could decay to the heavy E and its scalar 
superpartner E. This decay must compete with other decays into light fermions, however, 
and is likely be phase space suppressed. Whether it is observable depends on the details of 
the spectrum. Due to the approximate Z 2 symmetry, two Es must be produced at the end 
of any such cascade decay, which could give more lepton-rich signatures. 



C. Realization in the NMSSM 



Alternatively, for economy, one might try to embed the supersymmetric scenario in the 
NMSSM, adding only a pair of vectorlike leptons. The S field is already present to give the 
fj, term, and the lilghtest supersymmetric particle can provide the Dark Matter. Then we 
have: 

W N mssm = y£e R H d + ye R EH d + y s SEE + kS 3 + \SH u H d . (17) 

In this case X is identified with a mixture of higgsino and singlino. For the dark matter 
annihilation to be primarily into heavy leptons and not into W bosons, the dark matter 
must have a large singlino component. Achieving a cross section large enough to explain 



PAMELA requires the lightest pseudoscalar Higgs boson to be primarily singlet 67]. 
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VII. CONCLUSIONS 



The recent evidence of an excess of high energy positrons from the PAMELA experiment 
raises the tantalizing possibility that we may be seeing the first hints of a dark matter signal. 
If so, the hard lepton spectrum combined with the absence of any antiproton excess calls 
out for new thinking in terms of the dark matter annihilation processes. 

We have proposed a simple scenario which realizes the positron excess in PAMELA 
without exceeding the antiproton data. This required the introduction of a new heavy 
vectorlike lepton state that mixes with the lighter Standard Model leptons. We found that 
in addition to reproducing the PAMELA signal, this vectorlike lepton could be produced in 
Drell-Yan production at the LHC. Even if no missing energy signals are seen at the LHC, 
the observation of a vectorlike lepton, when combined with astrophysical observations might 
still illuminate the dark matter sector. 
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See, however, 6l| for a discussion of whether uncertainties in propagation parameters might 

accommodate annihilation to gauge bosons. 

We thank Tom Rizzo for a discussion of this point. 

A similar Lagrangian can be constructed for the situation where H is neutral (see, for instance, 



62fl). 

Here the boost refers specifically to the cross section into i/E rather than the total cross section. 
It should be noted that such a large cross section will necessitate some non-thermal mode of 
production of the dark matter to match the observed relic density. 
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